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ABSTRACT 


The  formulation  of  constitutive  relations  for  use  in  computer¬ 
ized  analyses  of  free-field  ground  shock  phenomena  is  based  primarily 
on  laboratory-determined  material  properties.  These  properties,  as 
described  by  stress-strain  relations ,  are  not  directly  determined  in 
the  laboratory,  but  are  derived  through  interpretation  of  load  and 
deformation  data  measured  by  the  experimenter.  Throughout  this  paper, 
one  laboratory  test,  the  triaxial  shear  test,  is  used  to  illustrate 
the  extent  of  interpretation  required  on  raw  data  and  the  influence 
of  this  interpretation  on  recommended  constitutive  properties.  Vari¬ 
ous  techniques  that  have  been  developed  to  obtain  stress-strain  data 
from  the  triaxial  test  are  reviewed  along  with  current  advances  in 
measurement  systems. 

Typical  raw  data  are  presented  and  calculations  of  axial,  lat¬ 
eral,  and  volumetric  strains  are  made  based  on  a  variety  of  empirical 
and  theoretical  approaches.  The  results  demonstrate  that  research 
and  development  efforts  are  still  required  in  the  area  of  material 
property  testing  in  order  to  establish  adequate  confidence  in  the 
formulation  of  constitutive  relations  for  ground  shock  calculations. 


PREFACE 


This  paper  was  prepared  for  presentation  at  the  Eric  H.  Wang 
Symposium  on  Protective  Structure  Technology  held  at  the  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico,  21-23  July 
1970.  The  subject  matter  presented  herein  was  primarily  intended 
for  those  persons  involved  in  the  field  of  ground  motion  prediction, 
but  not  necessarily  familiar  with  the  area  of  material  property 
determination . 

The  laboratory  equipment  and  techniques  described  in  this  re¬ 
port  were  developed  in  support  of  research  on  propagation  of  ground 
shock  through  soil  and  rock  being  conducted  by  personnel  of  the 
Soils  Division,  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES),  for  the  Defense  Atomic  Support  Agency  (DASA). 

This  report  was  prepared  and  presented  by  Mr.  J.  Q.  Ehrgott, 
Impulse  Loads  Section,  Soil  Dynamics  Branch,  Soils  Division,  WES. 
Helpful  comments  and  guidance  were  provided  by  Mr.  J.  G.  Jackson, 
Jr.,  Chief,  Impulse  Loads  Section.  Mr.  R.  W.  Cunny  was  Chief  of 
the  Soil  Dynamics  Branch  and  Mr.  James  P.  Sale  was  Chief  of  the 
Soils  Division.  Directors  of  the  WES  were  COL  Levi  A.  Brown,  CE, 
and  COL  Ernest  D.  Peixotto,  CE.  Technical  Director  was  Mr.  F.  R. 
Brown. 
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Area;  also,  calibration  error 

Original  area;  also,  fixed  area  of  ends  of  test  specimen 
Current  area  at  midheight  of  test  specimen 
Output  error 

Overall  accuracy  of  measurement  device 
Diameter  of  test  specimen 
Original  diameter  of  test  specimen 

Current  diameter  of  test  specimen  at  midheight  during 
hydrostatic  test 

Current  diameter  of  test  specimen  at  midheight  during 

shear  test 

Young's  modulus 

Shear  modulus 

Height  of  test  specimen 

Original  height  of  test  specimen 

Current  height  of  test  specimen  during  hydrostatic  test 

Current  height  of  test  specimen  during  shear  test 

Bulk  modulus 

Mean  normal  stress 

Operator  error 

Volume  of  test  specimen 


Current  volume  of  test  specimen  during  hydrostatic  test 
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Vq  Original  volume  of  test  specimen 

Corrected  original  volume  of  test  specimen 

Deformed  volume  of  test  specimen  at  end  of  hydrostatic  test 

Vg  Deformed  volume  of  test  specimen  during  shear  test 

AV/Vq  Volumetric  strain 

AV  Volume  change 

e  Axial  strain 

a 

e  Radial  strain 

r 

€q  Tangential  strain 

e  -  e  Deviator  strain 

a  r 

o  Axial  stress 

a 

or  Radial  stress;  also,  confining  pressure 

Og  Tangential  stress 

a  -  cr  Deviator  stress 

a  r 


CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 

British  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  units  as  follows. 


Multiply 


By 


To  Obtain 


mils  0.0254 

inches  2.54 

•feet  0 . 3048 

pounds  per  square  inch  0.070307 


millimeters 

centimeters 

meters 

kilograms  per  square 
centimeter 


pounds  per  cubic  foot 


16.01846 


kilograms  per  cubic  meter 


CHAPTER  1 
INTRODUCTION 


The  response  of  any  land-based,  structure  to  loadings  produced  by 
the  explosive  impact  resulting  from  either  conventional  high  explosive 
(HE)  or  nuclear  detonations  is  highly  dependent  upon  the  dynamic 
stress-strain  and  strength  characteristics  of  the  surrounding  and 
supporting  earth  material.  Current  wave-propagation  computer  codes 
used  in  the  prediction  of  free-field  stresses  and  motions  require 
stress-strain  constitutive  relations  for  the  earth  materials  at  the 
sites  of  interest.  Generally,  such  relations  are  derived  from  labo¬ 
ratory  tests  on  undrained  soil  and/or  rock  specimens  conducted  under 
a  variety  of  states  of  impulsive -type  stress  and  at  magnitudes  closely 
simulating  expected  field  levels.  In  cases  where  current  equipment 
limitations  only  allow  the  application  of  static  loadings,  extrapola¬ 
tion  of  test  results  must  be  made  to  reflect  dynamic  conditions. 

The  derivation  of  stress-strain  properties  from  laboratory  data, 
however,  is  subject  to  analysis  and  interpretation  by  the  experi¬ 
menter.  It  is  the  intent  of  this  paper  to  illustrate,  by  examples, 
the  influence  that  analyses  and  interpretations  have  on  stress- strain 
relations  obtained  from  raw  test  data.  It  is  realized  that  a  con¬ 
tinuing  effort  is  being  made  to  solve  and/or  improve  the  uncertainties 
affecting  laboratory  tests;  however,  the  techniques  presented  herein 
are  representative  of  the  approaches  currently  being  used  to  obtain 
stress-strain  relations  during  production- type  testing  programs. 

In  order  to  best  illustrate  how  the  raw  data  are  analyzed,  one 
laboratory  test,  the  triaxial  shear  test,  will  be  used  throughout 
this  paper;  similar  illust rat ions  can  be  obtained  from  any  other  lab¬ 
oratory  test  based  on  its  peculiar  limitations  and  problems.  Since 
many  of  those  involved  with  the  development  of  constitutive  relations 
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may  be  unfamiliar  with  seme  areas  of  laboratory  techniques,  Chapter  2 
of  this  paper  will  present  a  general  review  of  the  triaxial  test  and 
recent  advances  made  in  measurement  systems.  In  Chapter  3>  various 
methods  used  to  calculate  volumetric  strain,  deviator  stress,  and 
deviator  strain  from  raw  data  will  be  derived,  and  in  Chapter  4, 
stress-strain  relations  derived  using  the  various  calculation  methods 
and  raw  data  obtained  from  two  triaxial  tests  will  be  compared  and 
discussed. 


12 


CHAPTER  2 
THE  TR1AXIAL  TEST 


2.1  TEST  DESCRIPTION 

The  triaxial  test,  when  properly  instrumented  so  as  to  provide 
complete  load  and  defo  mation  data,  is  oie  method  for  determining 
soil  and  rock  stress-strain  behavior  and  shear  strength.  Because 
triaxial  test  equipment  is  designed  to  permit  separate  control  of 
both  lateral  and  axial  loadings,  the  test  can  provide  data  on  the 
fundamental  response  characteristics  of  soil  and  rock  under  a  wide 
variety  of  controlled  states  of  stress .  In  addition  to  stress  con¬ 
trol,  the  test  also  permits  control  of  loading  rate,  drainage  condi¬ 
tions,  and  specimen  size.  A  brief  history  of  the  development  of  the 
triaxial  test  is  contained  in  ASTM  Special  Technical  Publication 
Number  361  (Reference  l) ,  and  a  detailed  description  of  the  apparatus 
and  test  procedure  is  contained  in  Reference  2. 

Separate  pressure  control  systems  are  necessary  for  application 
of  the  several  possible  axial  and  lateral  stress  paths .  Application 
of  the  axial  lead  can  be  accomplished  by  any  number  of  methods  de¬ 
pending  on  whether  the  test  is  desired  to  be  stress- controlled  or 
strain-controlled.  When  the  specimen  is  to  be  loaded  to  failure,  a 
strain- control  method  should  be  employed;  when  the  behavior  of  the 
specimen  is  to  be  studied  at  less -than- failure  stress  levels,  a 
stress- control  method  is  preferable  because  of  theNregulation  required 
in  loading  increment.  Dead  loading  the  sample,  eitheK directly  or  by 
lever  systems,  is  probably  the  oldest  and  simplest  stress-control 
method.  Pneumatic  systems  that  apply  air  pressure  to  a  movable  piston 
can  be  used  to  develop  not  only  very  large  loads  (by  varying  the  pis¬ 
ton  area  ratio)  but  also  to  develop  rapid  loading  rates .  Hydraulic 
systems  employ  the  same  principle  but  are  best  suited  for  strain- 
control  testing.  The  use  of  motor-driven  gears  is  also  a  very 
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effective  method  of  applying  a  strain-controlled  load.  Vibratory 
loads  can  be  achieved  by  means  of  a  rotating  cam  system,  a  pneumatic 
system  employing  a  sequence  valve,  springs  with  a  deadweight  loading 
device,  and  electrically  controlled  devices.  Shock  loadings  can  be 
accomplished  by  dropping  weights ,  pneumatic  loading  systems ,  or  ex¬ 
plosive  charges. 

Confining  pressures  can  be  achieved  by  means  of  any  number  of 
systems  employing  air  compressors,  hydraulic  pumps,  bottled  gas,  or 
piston- type  multipliers.  At  lower  pressure  levels  (less  than  500 
psi1) ,  the  use  of  gas  as  a  chamber  fluid  provides  easy  control  of 
the  confining  pressure.  At  higher  pressures,  a  hydraulic  fluid  must 
be  used  for  reasons  of  safety;  however,  a  hydraulic  fluid  becomes 
difficult  to  control  as  the  specimen  changes  in  volume  and  a  more 
complex  supply  system  is  required. 

In  the  standard  triaxial  configuration,  the  specimen  is  placed 
in  a  cylindrical  chamber.  The  axial  load  is  applied  to  the  specimen 
by  a  piston,  which  enters  through  a  sealing  device  (piston  guide)  in 
the  top  of  the  triaxial  cell.  The  inside  of  the  chantoer  can  be  pres¬ 
surized  to  provide  a  lateral  loading  or  confining  pressure  on  the 
specimen.  Generally,  the  cylindrical  specimen  is  sealed  from  the 
chamber  fluid  by  a  thin  membrane,  and  the  specimen  is  sandwiched  be¬ 
tween  two  rigid  plates  that  provide  the  transfer  of  load  between  the 
piston  and  the  specimen. 

The  triaxial  specimen  has  a  cylindrical  shape  with  a  height-to- 
diameter  ratio  of  approximately  2:1.  Axial  load  is  applied  in  the 
vertical  direction  or  along  the  z  axis.  Lateral  loading  is  applied 
radially  around  the  specimen  or  along  the  r  axis.  In  the  general 
analysis  of  the  triaxial  specimen,  it  is  assumed  that  the  tangential 
stress  and  strain  are  equal  in  magnitude  to  the  radial  stress  and 

1  A  table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  10. 
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strain  (i.e.,  ct„  =  a  ,  e„  =  e  )  and  that  the  stress  and  strain  dis- 
trihution  throughout  the  specimen  is  uniform.  Although  investiga¬ 
tions  such  as  those  in  References  3  and  4  have  shewn  that  stress  and 
strain  distribution  within  triaxial  specimens  is  generally  not  uni¬ 
form,  the  uniformity  assumptions  are  currently  necessary  to  the  over¬ 
all  analysis  of  the  test  results. 

Since  the  specimen  is  completely  within  the  chanter  during  load¬ 
ing,  the  confining  pressure  not  only  acts  radially  on  the  speci¬ 

men,  but  also  vertically.  Hence,  the  axial  loading  on  the  specimen 
is  the  sum  of  the  axial  force  of  the  piston  p  and  the  axial  force 
exerted  by  the  confining  pressure .  The  total  axial  force  divided  by 
the  specimen  area  (perpendicular  to  direction  of  loading)  is  defined 

as  the  axial  stress  a  .  Soil  mechanists  define  the  difference  be- 

a 

tween  the  axial  si  ess  and  the  confining  pressure  as  the  deviator 
stress  (ct  -  a  ) . 

9,  I* 

Assuming  for  the  moment  that  the  test  specimen  is  elastic,  homo¬ 
geneous,  and  isotropic  and  deforms  as  a  cylinder.  Figure  2.1  illus¬ 
trates  some  of  the  states  of  stress  that  can  be  imposed  on  the  speci¬ 
men  along  with  typical  stress-strain  responses.  For  a  hydrostatic 

o  +  2ct 
&  r 

state  of  stress  where  ct  =  a  ,  the  mean  normal  stress  p  = - - - 

a  r  5  3 

increases  and  the  deviator  stress  (ct  -  ct  )  remains  zero.  The  spec- 

a  r 

imen  response  is  usually  plotted  as  a  pres sure -volumetric  strain 
curve,  whose  slope  is  the  bulk  modulus  k  .  In  a  shear  test,  the 
specimen  is  first  loaded  hydrostatically  to  some  level  of  p  .  In 
cases  where  the  confining  pressure  is  held  constant  while  the  devia¬ 
tor  stress  is  increased,  the  specimen  response  is  usually  plotted  as 
a  deviator  stress  versus  axial  strain  €  curve,  the  slope  of  which 
is  Young's  modulus  E  .  In  cases  where  the  confining  pressure  is 
decreased  while  the  axial  stress  is  increased  so  that  the  mean  normal 
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stress  is  held  constant,  the  response  is  usually  plotted  as  a  de¬ 
viator  stress  versus  deviator  strain  (c  re)  curve,  the  slope  of 
which  is  two  times  the  shear  modulus  G  .  In  addition  to  the  modu¬ 
lus  data,  yield  strength  can  also  he  obtained.  Several  shear  tests 
conducted  at  different  levels  of  p  provide  data  to  describe  a  yield 
envelope . 

Of  course,  many  other  variations  can  be  conducted,  such  as  the 
extension  test  in  which  the  confining  pressure  is  held  constant  and 
the  axial  stress  is  decreased.  Yield  strength  values  obtained  from 
such  tests  describe  a  lower  bound  yield  envelope. 

2.2  MEASUREMENT  SYSTEM 

At  the  heart  of  any  test  device  is  the  measurement  system  that  re¬ 
ports  the  specimen's  response.  In  the  triaxial  test,  as  in  most  tests 
designed  to  obtain  earth  material  properties,  there  is  no  one  universal 
measurement  system  because  of  the  broaa  range  of  loads  and  deformations 
encountered.  The  system  that  will  be  presented. is  unique  to  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  and  used  mainly  for 
soils  and  soft  rocks;  several  individual  units  are  used  for  various 
measurement  ranges. 

The  static  confining  pressure  felt  by  the  specimen  can  be  meas¬ 
ured  by  a  number  of  commercially  produced  items  such  as  pressure 
transducers  and  gages.  In  a  dynamic  test,  special  pressure  trans¬ 
ducers  are  mounted  within  the  confining  chamber. 

The  social  load  is  best  measured  directly  on  the  soil  specimen  to 
eliminate  the  influence  of  piston  friction,  especially  during  high- 
pressure  cyclic  tests.  The  internal  load  cell  must  be  precalibrated 
for  pressure  effects,  however.  In  a  dynamic  test,  two  load  cells  are 
used  (one  above  and  one  beneath  the  specimen)  to  provide  a  cross  ref¬ 
erence  as  well  as  to  monitor  wave  propagation  in  those  tests  in  which 


M 


the  specimen  migit  be  loaded  too  fast. 

Although  the  measurement  of  load,  and  pressure  is  fairly  straight¬ 
forward,  measurement  of  specimen  deformation  is  not  so  simple.  All 
measurements  on  undrained  specimens  must  be  made  external  to  the  speci¬ 
men,  and  the  devices  used  for  such  measurements  must  not  restrict  or 
otherwise  influence  the  natural  deformation  of  the  specimen.  Further, 
they  must  be  relatively  free  from  pressure  influences  and,  in  the  case 
of  dynamic  tests,  able  to  respond  to  the  speeds  involved.  Ihe  sim¬ 
plest  method  for  measuring  axial  deformation  is  to  measure  the  move¬ 
ment  of  the  axial  load  piston  outside  the  chamber.  However,  during 
the  hydrostatic  phase  of  a  triaxial  test,  a  problem  develops,  espe¬ 
cially  in  a  dynamic  test,  as  to  how  to  ensure  piston  contact  with  the 
specimen.  An  axial  deformation  method  recently  incorporated  into  the 
WES  dynamic  triaxial  test  device  employs  two  vertical  deformeters, 
consisting  of  commercially  produced  LVDT's  mounted  at  180  degrees  to 
each  other,  to  monitor  the  movement  of  the  top  cap.  Ihe  devices  al¬ 
low  for  free  deformation  of  the  specimen,  require  only  a  nominal 
amount  of  force  to  be  moved,  and  permit  in-place  calibration  and  pres¬ 
sure  checks.  An  external  piston  measurement  system  is  also  employed 
during  the  shear  phase  as  a  cross-check  of  systems,  since  the  appli¬ 
cation  of  deviator  stress  during  this  phase  ensures  contact  between 
the  piston  and  the  specimen. 

Ihe  lateral  deformations  of  the  specimen  are  very  difficult  to 
evaluate  properly.  Not  only  does  the  specimen  deform  laterally,  but 
also  downward,  so  that  points  on  its  surface  do  not  move  in  a  per¬ 
fectly  radial  plane.  Further,  the  specimen  does  not  necessarily  de- 
foim  uniformly  so  that  a  single  lateral  measurement  may  not  be  com¬ 
pletely  representative.  Coupled  with  these  problems  is  the  fact  that 
there  is  a  rubber  membrane  surrounding  the  specimen  and  its  contribu¬ 
tion  to  the  measured  movement  also  has  to  be  evaluated.  Hence,  in 
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order  to  make  any  measurement  of  strain  in  the  radial  direction,  some 
problems  have  to  be  neglected  during  the  actual  test  measurements  and 
then  accounted  for  during  evaluation  of  results. 

Of  the  many  devices  which  have  been  developed  for  measurement  of 
lateral  strain,  the  one  described  in  Reference  5  appears  to  be  the 
best  suited  for  static  tests.  The  device  measures  changes  in  diameter 
of  the  specimens  by  means  of  cantilever  springs,  which  are  instru¬ 
mented  with  electrical  resistance  strain  gages.  Hie  lateral  defor- 
meter  is  simple  in  design,  easy  to  install,  and  exerts  little  if  any 
force  on  the  specimen.  Several  of  these  devices  can  be  arranged  at 
different  elevation  levels  around  the  specimen  so  that  a  fairly  com¬ 
plete  deformation  profile  can  be  obtained.  Also,  and  most  important, 
the  device  allows  for  in-place  calibration,  including  that  for  pres¬ 
sure  effects.  Another  system,  utilizing  commercially  produced  LVDT's, 
is  best  employed  for  tests  in  which  the  confining  pressure  is  dynami¬ 
cally  applied.  This  device  has  physical  limitations,  such  as  range, 
and  it  is  easily  damaged  during  uncontrolled  specimen  rupture. 

Figure  2.2  shows  a  cutaway  schematic  of  the  WES  high-pressure 
triaxial  test  chamber  with  the  measurement  system  in  place  around  the 
specimen.  Hie  electrical  outputs  from  the  various  systems  go  through 
amplifiers  and  are  recorded  on  a  direct-writing,  light-beam 
oscillograph. 

2.3  MEASUREMENT  ERRORS 

For  a  given  measurement  device,  there  are  generally  three  sources 
of  error,  i.e.  calibration  error,  output  error,  and  operator  error. 
Calibration  error  (+A)  is  related  to  the  accuracy  of  the  device  used 
as  a  standard  and  the  associated  errors  that  occur  during  the  calibra¬ 
tion;  +2  percent  is  reasonable  in  a  product ion -type  program.  Output 
error  (+B)  is  a  function  of  the  measuring  device  and  system  including 
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such  items  as  electronics,  temperature  effects,  etc.  A  value  of  +1 
percent  is  reasonable;  however,  that  error  increases  to  about  +2  per¬ 
cent  when  operating  at  less  than  10  percent  of  full  range.  Operator 
error  (+R)  is  a  function  of  the  recorded  output  and  is  a  given  per¬ 
centage  of  the  range  of  current  interest.  For  an  oscillograph  rec¬ 
ord,  the  error  is  based  on  the  maximum  working  range  of  the  oscil¬ 
lograph  and  the  operator's  readability;  a  value  of  +0.5  percent  is 
considered  reasonable.  The  overall  accuracy  (+C)  of  the  measurement 
device  can  then  be  expressed  by  the  following  equation 


=  +  £(+B  x  full  electrical  output  value)" 


+  (+R  x  full  electrical  output  value)" 


Where: 


21  1/2 

+  (+A  X  measurement  of  interest)  J  (2.1) 


B  =  +1  to  2  percent 
R  =  jO.5  percent 
A  =  +2  percent 

For  example,  assume  that  a  test  is  to  be  conducted  during  which 
the  maximum  deformation  of  the  material  is  estimated  to  be  500  mils 
and  that  the  most  suitable  measurement  instrument  with  a  range  ex¬ 
ceeding  the  expected  maximum  deformation  has  a  full  linear  range  of 
1,000  mils.  Also  assume  that  the  recorded  output  of  the  unit  can  be 
increased  by  adjusting  the  amplifier  gain  so  that  the  full  electrical 
output  is  generated  by  a  500-mil  mechanical  input.  The  accuracy  of 
a  measured  500-mil  deflection  would  be 

T  X/2 

C  =  [(0.01  x  500)2  +  (0.005  X  500)2  +  (0.02  X  500)2] 


«  +11-1/2  mils  or  +2-l/3  percent 


If,  however,  the  actual  maximum  deformation  of  the  specimen  is  only 
50  mils  or  even  if  the  specimen  deformed  500  mils,  hut  data  are  re¬ 
quired  at  the  initial  part  of  the  test,  say  at  the  50-mil  deformation 
point,  then  the  accuracy  of  the  measured  50-mil  deformation  would  be 


C 


[( 


0.02  X  500)2  +  (0.005  X  500)2 


+  (0.02  X 


w  +10-1/3  mils  or  +21  percent 


It  can  be  noted  from  the  above  examples  that  care  must  be  taken 
to  select  measurement  instruments  that  provide  the  greatest  possible 
accuracy  over  the  range  of  interest.  The  selection  depends  on  two 
factors:  experience  in  estimating  material  response  behavior  and  a 
predetermined  measurement  objective  such  as  to  obtain  initial  stress- 
strain  data  or  to  obtain  data  at  the  maximum  or  yield  condition. 

2.4  TYPICAL  RESULTS 

The  recorded  load  and  deformation  time  histories  are  converted, 
through  a  set  of  analytical  assumptions,  to  stress  and  strain  his¬ 
tories  and  displayed  for  further  property  analyses  in  the  form  of 
stress-strain  plots.  Two  typical  WES  static  high-pressure  triaxial 
test  data  plates  are  shown  in  Figures  2.3  and  2.4;  Figure  2.3  illus¬ 
trates  results  for  a  constant  p-type  test  conducted  on  a  sandy  clay 
specimen  and  Figure  2.4  illustrates  results  from  a  constant  o^-type 
test  conducted  on  a  clayey  siltstone.  The  sandy  clay  specimen  (Fig¬ 
ure  2.3)  was  first  loaded  hydrostatically  to  1,000  psi  and  then 
loaded  to  failure  in  shear  by  increasing  the  deviator  stress  while 
the  mean  normal  stress  was  maintained  constant.  The  siltstone  speci¬ 
men  (Figure  2.4)  was  hydrostatically  loaded  to  1,500  psi  and  then 
loaded  to  failure  in  shear  by  increasing  the  deviator  stress  while 
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the  confining  pressure  was  held  constant.  The  response  of  each  speci 
men  to  the  hydrostatic  loading  is  seen  as  a  plot  of  pressure  versus 
volumetric  strain,  and  the  response  to  the  shear  loading  is  seen  as 
a  plot  of  deviator  stress  versus  deviator  strain.  Pertinent  informa¬ 
tion  regarding  the  composition  properties  of  the  specimens  is  also 
shown  on  the  figures. 
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CHAPTER  3 

INTERPRETATION  OF  RAW  DATA 


'  Between  collection  of  the  raw  test  data  and  the  production  of 
data  plots  such  as  those  shown  in  Figures  2.3  and  2.4,  much  inter¬ 
pretation  and  analysis  must  be  done  by  the  material  property  investi¬ 
gator.  This  analysis  of  the  raw  data  has  considerable  influence  on 
the  constitutive  properties  finally  selected  to  describe  the  mate¬ 
rial's  behavior.  The  stress-strain  plots  shown  in  the  above- 
mentioned  figures  were  based  on  a  certain  set  of  assumptions;  other 
stress-strain  plots  could  have  been  generated  from  the  same  set  of 
raw  data  using  different,  yet  perhaps  equally  valid,  assumptions. 

In  the  following  section,  some  of  the  various  methods  and  assumptions 
that  can  be  used  to  analyze  raw  data  will  be  described.  Then  the  data 
from  the  two  previously  presented  triaxial  tests,  one  on  a  relatively 
soft  sandy  clay  (Figure  2.3)  and  one  on  a  relatively  stiff  clayey 
siltstone  (Figure  2.4),  will  be  used  to  illustrate  the  variations  in 
moduli  that  can  be  obtained  using  the  different  methods  of  analysis. 

3.1  DETERMINATION  OF  VOLUMETRIC  STRAIN 

Since  deformation  measurements  are  made  only  at  a  few  specific 
locations  on  the  specimen  surface,  it  is  necessary  to  assume  a  com¬ 
plete  specimen  deformation  pattern  in  order  to  describe  the  overall 
stress-strain  response  of  the  specimen  to  imposed  loadings.  First, 
consider  some  of  the  general  deformed  shapes  noted  during  application 
of  hydrostatic  loadings.  Figure  3 »1  shews  the  two  most  common  shapes, 
double  cone  and  cylindrical.  Specimens  deforming  as  a  double  cone 
(IA,  IB,  IC)  undergo  maximum  lateral  deformation  at  the  center  while 
the  ends  undergo  little  or  no  lateral  deformation.  The  deformation 
restraint  at  the  ends  is  attributed  to  end-cap  friction  (References 
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6  and  7)  •  Shear  forces  at  both  ends  caused  by  the  top  cap  and  base 
pedestal  prevent  free  movement  of  the  specimen;  however,  the  influ¬ 
ence  of  end  caps  is  not  the  same  on  all  types  of  material.  Figures 
II-A,  II -B,  and  II-C  of  Figure  3»1  show  the  second  type  or  cylin¬ 
drical  shape.  The  specimens  deform  radially  in  a  uniform  or  nearly 
uniform  manner.  The  effect  of  the  end  caps  on  the  specimen  is  not  as 
apparent  as  in  the  case  of  the  double -cone  shape.  There  is  little 
factual  data  regarding  the  distribution  of  axial  deformation  through¬ 
out  the  specimens  in  either  of  the  two  cases,  although  it  is  surmised 
that  the  distribution  is  more  nearly  uniform  in  the  cylindrical  shape. 

Figures  3*2  and  3*3  show  two  separate  recompacted  clayey  silt 
specimens  after  being  subjected  to  confining  pressures  of  500  and 
5,000  psi,  respectively.  Although  some  rebounding  occurred,  the 
general  deformed  shape  achieved  after  completion  of  loading  remained 
upon  removal  of  the  pressure.  The  shape  is  basically  double  cone 
with  some  deformation  of  the  ends  as  typified  in  Figure  I-C  of  Fig¬ 
ure  3*1«  Figures  3.4  and  3.5  show  deformed  specimens  of  silty  clay 
with  rock  fragments  after  hydrostatic  loading  to  pressures  of  1,000 
and  5 >000  psi,  respectively.  Note  that  these  specimens  deformed  as 
fairly  uniform  cylinders  as  typified  in  Figure  II-A  of  Figure  3*1* 

All  of  the  above  specimens  were  tested  in  the  same  manner;  thus  the 
differences  in  deformed  shape  can  only  be  attributed  to  the  physical 
properties  of  the  specimen  material.  Each  type  material,  therefore, 
should  be  considered  separately  in  analyses. 

To  illustrate  the  various  approaches  that  can  be  used  to  calcu¬ 
late  volumetric  strain,  consider  only  one  of  the  deformed  shapes  as 
shown  in  Figure  3*1*  In  this  case  the  specimen  starts  from  an  ini¬ 
tial  shape  of  a  -uniform  right  circular  cylinder  whose  volume  V 


Vo  -  f  <V2  <«o> 


(3-1) 


Where: 


D  =  original  diameter 
o 

H  =  original  height 
o 

The  specimen  then  deforms  as  a  double  cone  with  fixed  ends,  i.e., 
shape  I-B.  Throughout  the  example,  compression  will  be  considered  as 
positive  and  all  strains  will  be  calculated  in  Lagrangian  notation, 
i.e.  (change  in  dimension)/( original  dimension).  Raw  data  measure¬ 
ments  include  pressure  p  ,  total  height  change  AH  ,  and  total  diam¬ 
eter  change  at  midheight  of  the  specimen. 

3.1.1  Method  V-l.  In  this  method,  the  deformed  shape  is  approx¬ 
imated  by  straight  lines  and  hence,  two  truncated  cones  placed  end  to 
end  as  shown  in  Figure  3*6.  The  volume  Vc  of  the  specimen  at  any 
given  time,  i.e.  the  current  volume  is 

(3-2) 


Where: 


V  =  initial  volume  =  H  A 
o  o  o 

Hn  =  current  height  =  H  -  AH  _2 
1  O  TTD^ 

=  current  area  at  midheight  =  -jp 
=  current  diameter  at  midheight  =  D  -  AD 

A  =  fixed  area  of  ends  =  — 


The  volumetric  strain  —  can  be  found  from 

o 

AW  V  -  V 
AV  o  c 


?  DX  -  f  (£)  (Di +  Do +  Vo) 


?d2h 
4  o  c 


(3.3) 
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Substitutions  can,  of  course,  be  made  to  allow  comparison  of 
volumetric  strain  in  terms  of  axial  and  radial  strains,  i.e. 


€  =  axial  strain  =  ^ 

a,  n 

O 

(3.4) 

e  =  radial  strain  =  — 
r  D 

0 

• 

(3.5) 

H,  =  H  -  e  H 

J-  0  a  0 

(3.6) 

D,  =  D  -  e  D 

J.  0  r  0 

(3.7) 

and  the  volumetric  strain  for  Method  V-l  becomes 

M  -  ,BA  ‘  (Ho  "  Vo’  k  k  -  €rD0)2  +  £  +  (»„  - 

v>„>  <vl 

€2 

=  £a+er-  Va+f<e.-« 

(3.8) 

This  is  the  method  currently  used  to  calculate  the  volumetric 
strain  as  presented  in  Figures  2.3  and  2.4,  and  will  be  subsequently 
referred  to  as  the  standard  method.  This  does  not  imply  that  it  is 
necessarily  the  correct  method,  but  rather  only  a  standard  used  in 
this  paper  for  comparison  purposes. 

3.I.2  Method  V-2.  In  this  method,  the  specimen  is  assumed  to 
always  have  an  average  diameter  of  (Dq  +  D^/2  as  indicated  in  Fig¬ 
ure  3.7.  The  volume  Vq  of  the  specimen  at  any  given  time  is 


Vc  *  Vo  -  AV  ■?(»„-  f  ^  (3.9) 
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where  AD^  =  current  total  change  in  diameter  at  midheight;  other 
notations  are  as  before.  Hie  volumetric  strain  can  be  found  from 


AV 


2  . 

g  fo,  ~  g  (Bo  ~  I  *>l)  h 


V°  D^o 

Keeping  the  same  definition  of  axial  and  radial  strain,  the  volu¬ 
metric  strain  for  Method  V-2  is  expressed  as 


(3.10) 


—■=  e.  +  e  -  ee  +  -jp  \e  -  l)  (3*ll) 

V  a  r  a  r  4  a 
o 

3.1.3  Method  V-3«  In  this  method,  the  midheight  lateral  de¬ 
formation  is  assumed  to  be  representative  of  the  entire  specimen  as 
indicated  in  Figure  3*8*  Hie  current  volume,  therefore,  is 


vc  -  To  -  Av  =  g  Dfr. 
Hie  volumetric  strain  can  be  found  from 


(3-12) 


(3-13) 


and  in  terms  of  axial  and  radial  strains,  the  volumetric  strain  for 
Method  V-3  is 


A l 

V 

o 


O 

e  +  2e  -  2e  e  +  e  (e  -  l) 
a  r  r  a  r  a 


(3.14) 
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3.1.4  Method  V-4.  In  Method  V-3,  the  current  volume  of  the  spec¬ 
imen  was  assumed  to  be  a  right  circular  cylinder  whose  diameter  was 
equal  to  the  midheight  diameter,  and  the  volume  outside  that  cylinder 
caused  by  end  cap  friction,  i.e.  dead  zone, was  neglected.  In  this 
method,  the  current  volume  is  the  same,  but  the  original  volume  of 
the  specimen  is  corrected  to  reflect  a  loss  in  volume  due  to  the  pre¬ 
viously  neglected  dead  zone  as  shown  in  Figure  3-9*  The  corrected 
original  volume  of  the  specimen  at  any  given  time  is 


ro  -  5  Dk  -  [t  (?)  (D? +  Do +  Vo)  -  5  DX] 


Volume  of  Dead  Zone 


The  current  volume  of  the  specimen  is 


(3-15) 


vc  "  ?  <V2 


(3-16) 


Ip  this  method,  volumetric  strain  is  defined  as  the  change  in  volume 
divided  by  the  corrected  volume,  or 


V1  -  V 
AV  o  c 


(3.17) 


After  similar  substitutions  as  presented  in  the  previous  methods. 
Method  V-4  becomes 


e  +  e  -  e  e  -  i  e2  (l  -  e  ) 
a  r  a  r  3  r  a 

1  -  e  +  e  e  +  f  e2  (l  -  e  ) 
r  r  a  3  r  a 


(3.18) 


3.1.5  Method  V-5.  In  this  method,  it  is  assumed  that  only  an 
axial  deformation  measurement  is  available  and  also  that  all  strains 
are  equal  during  hydrostatic  loading  in  accordance  with  elastic  theory. 


Therefore 


f'3£a 

o 

3.1.6  Method  V-6.  This  method  is  similar  to  Method  V-5  except 
that  it  is  assumed  that  only  a  radial  deformation  measurement  is  avail¬ 
able.  Therefore 


&L  -  «  (3.20) 

V  “  ~3  r 
o 

3.1.7  Summary.  Other  calculation  methods  could  be  developed 
although  each  of  the  above  six  methods  has  been  based  to  some  extent 
on  observed  phenomena  and  elastic  theory.  In  Method  V-l,  the  volu¬ 
metric  strain  was  based  on  the  geometric  shape  of  the  specimen. 

Method  V-2  was  based  on  an  average  deformation  assuming  that  the  ra¬ 
dial  deformation  at  the  center  is  the  extreme  for  the  specimen. 

Method  V-3  considered  only  the  center  deformation  as  being  representa¬ 
tive  and  neglected  distortions  due  to  end  restraint.  Method  V-4  con¬ 
sidered  the  specimen  to  be  divided  into  two  zones,  a  center  cylindri¬ 
cal  zone  and  a  surrounding  dead  zone.  Methods  V-5  and  V-6  are  based 
entirely  on  theory  of  elasticity  and  the  assumption  that  only  one  of 
the  measured  deformations  is  valid.  The  consistent  Lagrangian  defi¬ 
nition  given  throughout  the  above  calculations  to  axial  and  radial 
strains  does  not  imply  that  they  are  true  strain  values  nor  that 
these  strains  are  distributed  uniformly  within  the  specimen.  The 
purpose  of  the  notation  was  only  to  permit  easy  comparison  c  the 
different  equations  of  volumetric  strain.  A  summary  of  the  methods 
is  presented  in  Figure  3*10. 


3.2  DETERMINATION  OF  DEVIATOR  STRESS  AND  STRAIN 

When  the  shear  phase  of  the  test  is  initiated,  the  specimen  has 
already  deformed  due  to  the  prior  hydrostatic  phase  to  one  of  the 
various  shapes  discussed  in  the  previous  section.  From  this  initial 
condition,  the  specimen  deforms  downward  axially  and  outward  radially 
under  the  influence  of  increased  deviator  stress  (ct&  -  cr^) .  Shape  is 
again  important  in  quantitatively  determining  the  strain  and  stress; 
however,  it  is  more  difficult  to  properly  evaluate  shape  during  the 
shear  phase  due  to  the  formation  of  a  Complicated  series  of  shear 
zones .  ! 

Various  shapes  observed  after  shear  failure  of  both  soils  and 
soft  rocks  are  depicted  in  Figure  3*U*  Shear  Types  IA,  IB,  and  IC 
represent  ductile-type  failures  with  a  bulge  as  the  predominant  fea¬ 
ture.  Shear  Types  IIA,  I  IB,  and  IIC  typify  more  brittle  materials 
with  a  shear  plane  as  the  predominant  feature. 

Type  IA  shows  the  typical  bulge-type  failure  with  the  influence 
of  end-cap  friction  preventing  deformation  of  the  specimen  at  the 
ends.  The  failure  shown  in  IB  indicates  the  uniform  deformation  ex¬ 
pected  of  a  ductile-type  material  free  of  end  restraint.  Type  IC 
shows  a  emiductile  material  in  which  a  bulging-type  deformation  dom¬ 
inates,  but  cracking  is  noted,  probably  due  to  end-cap  friction  caus¬ 
ing  the  formation  of  cone-shaped  zones  at  either  end  of  the  specimen. 

Type  IIA  is  the  classical  type  failure  with  a  definite  shear 
plane  formed  within  the  specimen;  such  a  failure  for  a  sandstone 
specimen  sheared  during  a  dynamic  constant  p-type  test  is  shown  in 
Figure  3*12.  The  failure  typified  by  I IB  indicates  the  influence  of 
end-cap  friction  with  a  distinct  cone  or  wedge  formed  at  either  or 
both  ends  of  the  specimen.  The  wedge-type  failure  is  usually  accom¬ 
panied  by  vertical  fracture  zones  along  the  axis  of  the  specimen. 

Type  IIC  illustrates  a  composite  type  failure  where  the  specimen, 
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composed  of  two  types  of  material,  fails  in  both  materials  or  in  just 
the  weaker  material. 

In  Figure  3*3 5  the  double-cone  shape  was  shown  as  formed  after 
application  of  a  5,000-psi  hydrostatic  loading  to  a  recontacted  clayey 
silt  specimen.  Figure  3*13  shows  another  specimen  of  the  same  mate¬ 
rial  after  a  small  increment  of  deviator  stress  has  been  applied  sub¬ 
sequent  to  a  5,000-psi  hydrostatic  loading.  Note  the  small  outward 
bulging  just  starting  to  occur  at  the  center  of  the  specimen  while 
the  ends  sure  relatively  unaffected. 

The  next  series  of  photographs  (Figures  3*14  through  3«17)  il¬ 
lustrates  development  of  a  Type  IA  failure.  Several  essentially 
identical  specimens  of  modeling  clay  were  first  subjected  to  a  con¬ 
fining  pressure  of  5j000  psi  and  then  to  various  increments  of  axial 
deviator  stress  up  to  the  peak  yield  strength.  Figure  3*14  shows  an 
original,  undeformed  specimen.  Figure  3-15  shows  a  specimen  after  a 
small  increment  of  axial  deviator  stress  has  been  applied  and  removed; 
the  center  bulges  outward  while  the  ends  remain  unaffected.  The 
photograph  in  Figure  3*16  was  taken  after  a  greater  application  of 
deviator  stress;  it  shows  that  the  center  bulge  area  has  increased 
while  the  ends  have  just  started  to  deform.  After  peak  yield  strength 
has  been  reached,  the  specimen  assumed  the  shape  shown  in  Figure  3«17- 
Apparent  dead  zones  occurred  at  either  end  of  the  specimen  along  with 
noticeable  deformations  at  the  quarter  heights  of  the  specimen.  The 
nonuniform  distribution  o.’  stress  and  strain  within  the  specimen  is 
obvious.  The  sketches  in  Figure  3*18  illustrate  the  various  deformed 
shapes  which  develop  during  a  complete  triaxial  shear  test. 

As  with  the  case  of  the  hydrostatic  test,  the  deformed  shape 
can  be  used  as  a  guide  to  calculate  the  average  stress  and  strain 
within  the  specimen.  To  illustrate  the  various  methods  that  can  be 
used  to  calculate  deviator  stress  and  strain,  it  will  be  assumed  that 
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measurements  of  deviator  load,  confining  pressure,  axial  deformation, 
and  center  radial  deformation  are  available  and  that  the  specimen  ex¬ 
hibited  a  double-cone  shape  (Type  IB)  during  hydrostatic  loading  and 
a  bulge-type  failure  (Type  IA)  when  sheared.  Strains  will  be  defined 
in  terms  of  measured  deformations  for  purposes  of  comparison. 


«  .ffe 

a  H 

o 

Where: 


AHg  =  change  in  height  measured  during  shear  phase 


Where: 


(3.21) 


(3-22) 


ADg  =  change  in  midheight  diameter  measured  during  shear  phase 

H  and  D  are  the  original  prehydrostatic  phase  height  and  diam- 
o  o 

eter,  respectively,  so  that  the  total  strain  from  the  start  of  the 
hydrostatic  phase  to  some  point  of  interest  during  the  shear  phase 
can  be  found  by  algebraic  addition  of  the  strains  from  both  phases. 

Ihe  radial  strain  during  the  shear  phase  is  negative  since  it  will 
be  an  outward  movement. 

3.2.1  Method  S-l.  In  this  method,  it  is  assumed  that  the  axial 
strain  and  the  radial  strain  as  calculated  from  the  vertical  displace¬ 
ment  of  the  top  cap  and  the  midheight  radial  deflection  completely 
represent  the  actual  strains  occurring  within  the  specimen.  Figure 
3.19  shows  the  assumed  deformed  shape  at  the  start  of  the  shear  test 
and  during  the  shear  test.  The  deviator  strain  becomes 


and  the  deviator  stress  is 


(3.23) 
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Where: 


=  current  diameter  during  shear  test 
This  method  was  used  to  calculate  the  deviator  stress  and  strain  re¬ 
sults  shown  in  Figures  2.3  and  2.4  and  will  subsequently  be  referred 
to  in  this  paper  as  the  "standard"  method  for  comparison  purposes. 

3.2.2  Method  S-2.  For  this  method,  assume  that  the  axial  de¬ 
formation  is  representative,  but  that  the  center  radial  deformation 
is  a  maximum  value  and  therefore  not  representative  of  the  entire 
specimen.  The  radial  deformation  should  therefore  be  weighted,  and 
for  simplicity,  an  average  value  of  one-half  is  used  as  shown  in  Fig¬ 
ure  3*20.  Hie  deviator  strain  becomes 


(3-25) 


and  the  deviator  stress  is 


(3.26) 


Where: 

D^  =  Deformed  diameter  at  the  end  of  the  hydrostatic  phase 
3.2.3  Method  S-3»  In  this  method,  it  is  assumed  that  only  the 
axial  measurement  is  available  and  that  the  ratio  of  radial  strain  to 
axial  strain  is  equal  to  -0.5  in  accordance  with  elastic  theory  for 
pure  shear.  The  deviator  strain  is  therefore 


e 

a 


(3.27) 
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and  the  deviator  stress  is 
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3.2.4  Method  S-4.  This  method  is  derived  based  on  the  assump¬ 
tions  that  only  the  measured  axial  deformation  is  valid  and  that  no 
volume  change  occurs  during  the  shear  phase.  This  method  is  commonly 
used  in  soil  mechanics  where  only  axial  deflection  is  measured.  Hie 
specimen  is  assumed  to  undergo  cylindrical  deformation  as  shown  in 
Figure  3*21  with  no  change  in  volume.  Because  it  is  assumed  that 
there  is  no  volume  change  during  shear,  the  current  diameter  during 
the  shear  test  can  be  expressed  as 


D 


2 


(3.29) 


and  the  radial  strain  can  be  calculated  as 


6 

r 


Where: 

=  deformed  volume  at  the  end  of  the  hydrostatic  phase 
Hg  =  current  height  during  shear  phase 
The  deviator  strain  is  therefore 
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(3.31) 
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and  the  deviator  stress  is  found  from 


a 

a 


(3.32) 


3.2.5  Method  S-5«  In  this  method,  as  in  Method  S-4,  it  is  as¬ 
sumed  that  there  is  no  volume  change,  but  in  this  case  the  axial  de¬ 
flection  is  assumed  to  be  invalid  so  that  only  the  radial  deformation 
is  available  for  calculation  purposes.  The  volume  of  the  specimen. at 
the  end  of  the  hydrostatic  test  V.^  and  the  current  diameter  during 
the  shear  phase  D ^  may  be  used  to  calculate  an  assumed  current 
height  during  the  shear  phase  HJ^  for  a  right  circular  cylindrical 
specimen,  i.e. 
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The  axial  and  radial  strains  may  be  expressed  as 
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(3.35) 


The  deviator  strain  is  therefore 


(3.36) 
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and.  the  deviator  stress  is 
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3 .2.6  Method  S-6.  For  this  method,  it  is  assumed  that  the  axial 
strain  as  calculated  from  the  measured  change  in  specimen  height  is 
not  representative  of  the  axial  strains  within  the  specimen  and  that 
the  change  in  height  measurement  should  he  corrected  by  some  empiri¬ 
cal  factor.  An  arbitrary  factor  of  2  is  used  for  illustration  pur¬ 
poses.  The  deviator  strain  in  this  case  is 
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and  the  deviator  stress  becomes 
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(3.39) 


3.2.7  Summary.  The  above  six  methods  represent  various  proce¬ 
dures  for  calculating  stress  and  strain  in  soil  specimens  during  the 
shear  phase  of  triaxial  tests.  In  Method  S-l,  stress  and  strain  were 
based  on  actual  measurements  of  applied  loads,  displacement  of  the  top 
cap,  and  midheight  diameter  changes.  Method  S-2  was  based  on  a  repre¬ 
sentative  diameter  as  the  average  of  a  fixed  end  diameter  and  the  meas¬ 
ured  diameter  at  the  center  of  the  specimen.  Method  S-3  is  based  on 
an  assumed  strain  ratio  of  -0.5  for  elastic  pure  shear  and  the  assump¬ 
tion  that  only  the  axial  deformation  measurements  are  valid.  For  both 
Methods  S-4  and  S-5,  it  was  assumed  that  no  volume  change  occurs  dur¬ 
ing  the  shear  phase.  Method  S-4  required  the  use  of  only  the  axial 
measurement  for  calculations;  Method  S-5  required  only  the  use  of  the 
radial  measurement.  Method  S-6  was  an  attenqpt  to  correct  the  axial 
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deformation  by  means  of  an  arbitrary  empirical  factor  based  on  limited 
research  conducted  on  the  distribution  of  axial  deformation  of  undis¬ 
turbed  specimens  during  the  shear  phase.  In  general,  the  greatest 
axial  strain  seems  to  occur  near  the  midheight  of  the  specimen  as  il¬ 
lustrated  in  Figure  3.22.  All  six  methods  discussed  in  this  section 
are  summarized  in  Figure  3*23* 
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Figure  3.1  Deformed  shapes  of  specimens  during  hydrostatic  loading. 
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ORIGINAL  DEFORMED 


Ca  =  CTr  =  500  PSI 

- O 

Figure  3.2  Triaxial  specimen  of  recompacted  clayey  silt  after  being 
subjected  to  500-psi  hydrostatic  pressure,  shape  IC. 


cra=c^  =  5000  psi 


Figure  33  Triaxial  specimen  of  recompacted  clayey  silt  after  beii 
su  jected  to  5,000-psi  hydrostatic  pressure,  shape  IC. 
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.  DEFORMED  SHAPE  b.  ASSUMED  DEFORMED  c.  DOUBLE  CONE 

SHAPE 

Figure  3*6  Cross  section  of  shapes  considered  in  Method  V-l. 


Figure  3*7  Cross  section  of 
assumed  shape  used  in  Method 

V-2. 


Figure  3*8  Cross  sec¬ 
tion  of  assumed  shape 
used  in  Method  V-3. 


Figure  3*9  Assumed  shape 
used  in  Method  V-4 . 
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Figure  3-H  Shapes  of  failed  specimens  after  shear  test. 


b .  Specimen  separated  to  show  shear  plane . 

Figure  3.12  Triaxial  specimen  of  sandstone  after  shear  failure  during 
a  dynamic  constant  p-type  test,  shape  IIA. 
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Figure  3.13  Triaxial  specimen  of  recompacted  clayey  silt  after 
being  subjected  to  a  small  deviator  stress  while  maintaining  a 
confining  pressure  of  5,000  psi. 
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ORIGINAL  SHAPE 


Figure  3.14  Triaxial  specimen  of  modeling 
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Figure  3-15  Triaxial  specimen  of  modeling  clay  after  ap 
plication  of  a  small  deviator  stress  while  a  confining 
pressure  of  5,000  psi  is  maintained. 
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Figure  3-16  Tri axial  specimen  of  modeling  clay  after  ap¬ 
plication  of  a  larger  deviator  stress  vhile  a  confining 
pressure  of  5 >000  psi  is  maintained. 
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Figure  3*17  Triaxial  specimen  of  modeling  clay  after  ap¬ 
plication  of  large  (postyield)  deviator  stress  while  a 
confining  pressure  of  53000  psi  is  maintained. 
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Figure  3.18  Specimen  deformation  during  shear  test. 
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Figure  3.19  Cross  section  of  specimen  showing  assumed 
deformed  shapes  considered  in  Method  S-l. 


Figure  3*20  Cross  sec¬ 
tion  of  assumed  deformed 
shape  considered  in 
Method  S-2. 


Figure  3*21  Cross  section 
of  assumed  deformed  shape 
coisidered  in  Method  S-4. 
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o.  DEFORMED  SHAPE 
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b.  DISTRIBUTION  OF  AXIAL  STRAIN 


Figure  3.22  Distribution  of  axial  strain  along  length  of  a  specimen, 
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CHAPTER  4 

DISCUSSION  OF  INTERPRETATION  METHODS 


As  described  in  Chapter  3,  various  methods  can  be  used  to  calcu¬ 
late  stress  and  strain  from  load  and  deformation  measurements  obtained 
during  triaxial  tests.  The  differences  between  each  method  can  prob¬ 
ably  best  be  illustrated  by  stress-strain  curves  calculated  from  raw 
test  data  obtained  from  two  extremely  different  undisturbed  specimens, 
i.e.  one  quite  soft  and  one  very  stiff.  The  soft  specimen,  the  sandy 
clay  previously  described  in  Figure  2.3,  was  subjected  to  a  hydrostatic 
loading  of  1,000  psi  and  then  carried  to  failure  in  shear  while  the 
mean  normal  stress  remained  constant.  The  stiff  specimen,  the  silt- 
stone  previously  described  in  Figure  2.4,  was  subjected  to  a  hydro¬ 
static  loading  of  1,500  psi  and  then  carried  to  failure  in  shear  while 
the  confining  pressure  was  held  constant.  Measurements  made  for  both 
tests  included  confining  pressure,  axial  load,  change  in  height,  and 
change  in  midheight  diameter.  A  double- cone  deformed  shape  best  ap¬ 
proximated  the  response  of  both  specimens  during  the  hydrostatic  test, 
followed  by  a  bulging  shape  during  the  shear  test.  Both  specimens 
formed  shear  planes  at  failure.  As  previously  mentioned,  Method  V-l 
for  the  hydrostatic  test  and  Method  S-l  for  the  shear  test  will  be 
considered  as  the  standard  analysis  methods  for  purposes  of  comparison. 

4.1  COMPARISON  FOR  HYDROSTATIC  TESTS 

Pressure-volumetric  strain  plots  for  the  sandy  clay  specimen  dur¬ 
ing  the  hydrostatic  loading  test  are  shown  in  Figure  4.1  for  each  of 
the  six  methods  previously  described  for  calculation  of  voluhvetric 
strain.  Values  of  volumetric  strain  and  approximate  bulk  moduli  K  , 
calculated  for  both  the  0-  to  100-psi  and  the  600-  to  1,000-psi  pres¬ 
sure  ranges,  are  also  included  in  Figure  4.1.  Method  V-2  yielded 
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essentially  the  same  strain  values  as  the  standard  Method  V-l,  because 
of  the  equality  of  the  two  methods  except  in  small-order  terms . 

Method  V-4  resulted  in  only  slightly  greater  volumetric  strain;  i.e. 
the  difference  between  this  method  and  the  standard  is  unnoticeabls 
on  the  plots  for  volumetric  strains  less  than  5  percent,  which  is 
reasonable  to  expect  since  the  dead  zone  considered  in  the  calculation 
was  very  small  compared  to  the  large  original  volume.  Methods  V-6, 
V-3,  and  V-5  gave  substantially  larger  strains  than  did  the  standard 
method.  At  lower  pressures  (less  than  100  psi) ,  these  three  methods 
are  in  close  agreement  with  strains  exceeding  the  standard  by  approxi¬ 
mately  50  percent.  At  higher  pressures,  i.e.  1,000  psi,  Method  V-5 
gave  the  largest  strain,  exceeding  the  standard  by  approximately  65 
percent.  Methods  V-6  and  V-3  showed  fairly  close  agreement,  exceeding 
the  standard  strain  at  ,1,000  psi  by  about  45  percent.  It  should  be 
noted  that  since  the  deformation  results  are  material-response  depen¬ 
dent,  the  same  percentage  difference  may  not  hold  for  different  ma¬ 
terials.  In  fact,  similar  calculations  by  WES  using  other  test  data 
resulted  in  Method  V-6  giving  the  largest  volumetric  strain  values; 
however.  Methods  V-l,  V-2,  and  V-4  will  probably  always  give  smaller 
volumetric  strains  than  Methods  V-3,  V-5,  and  V-6. 

Comparison  of  bulk  moduli  values  for  the  pressure  range  of  0  to 
100  psi  shows  that  the  maximum  bulk  modulus  was  calculated  from  Method 
V-l.  The  lowest  modulus,  calculated  from  Method  V-6,  was  35  percent 
less  than  the  standard.  Values  of  bulk  moduli  for  the  high  pressure 
range,  600  to  1,000  psi,  varied  from  57  percent  less  than  to  21  per¬ 
cent  greater  than  the  standard.  The  lowest  modulus  was  calculated  by 
Method  V-5  and  the  greatest  modulus  by  Method  V-6.  Again  because  of 
the  dependence  of  the  deformations  on  material- response  behavior, 
there  are  no  fixed  relations  between  methods . 

The  same  basic  trends  observed  for  the  relatively  soft  sandy 
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clay  specimen  are  apparent  in  the  pressure-volumetric  strain  plots 
for  the  relatively  stiff  siltstone  specimen  as  shown  in  Figure  4.2. 
Methods  V-l,  V-2,  and  V--4  gave  essentially  the  same:  volumetric  strain 
results.  Again  Methods  ''-3,  V-5,  and  V-6  gave  substantially  larger 
strains  ranging  from  40  percent  (Method  V-6)  to  63  percent  (Method 
V-5)  greater  than  the  standard  at  a  pressure  of  1,500  psi.  At  a  pres¬ 
sure  of  20G  psi,  the  volumetric  strain  calculated  by  Method  V-5  was 
110  percent  greater  than  that  calculated  by  the  standard  method.  Note 
also  that,  while  Method  V-6  gave  larger  strains  at  high  pressure,  it 
gave  lower  strain  values  than  the  standard  method  in  the  low  pressure 
range . 

The  bulk  moduli  values  varied  from  a  low  of  42  percent  less  than 
that  calculated  from  the  standard  at  low  pressures  (200  to  600  psi)  ',0 
45  percent  less  thai  the  standard  at  high  pressures  (1,200  to  1,500 
psi) .  Method  V-5  gave  the  lowest  modulus  for  the  low  pressure  range 
while  Method  V-6  gave  the  lowest  value  for  the  high  pressure  range. 

At  the  low  pressure  range,  the  bulk  modulus  as  calculated  from  Method 
V-6  compared  favorably  with  the  standard. 

From  the  two  above  examples,  several  rather  general  observations 
can  be  made  regarding  the  comparison  of  methods  to  calculate  volu¬ 
metric  strain.  First,  the  results  indicate  that  the  volumetric 
strains  are  different  at  any  given  pressure  when  calculated  by  Meth¬ 
ods  V-5  and  V-6.  Both  methods  are  based  on  the  theory  of  elasticity, 
which  assumes  for  hydrostatic  compression  that  axial  and  radial 
strains  are  equal;  one  method  used  only  axial  strains  and  the  other 
only  radial  strains  to  calculate  volumetric  strains.  Obviously  the 
axial  and  radial  strains  as  derived  from  the  axial  and  radial  defor¬ 
mations  of  the  specimens  were  not  equal  nor  were  the  specimens 
elastic. 

Far  from  being  elastic  solids ,  soil  specimens  are  in  fact 
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multiphase  mixtures  of  air,  water,  and  mineral  particles.  Thus,  under 
hydrostatic  loading,  overall  specimen  compressibility  is  governed 
primarily  by  the  compressibility  of  the  air  as  it  is  driven  into 
solution.  After  full  saturation  has  been  achieved,  specimen  compress¬ 
ibility  is  then  governed  only  by  the  compressibility  of  the  water 
and  solids,  which  should  result  in  noticeable  stiffening  in  pressure 
versus  volumetric  strain  behavior.  The  volumetric  strain  at  which 
this  stiffening  occurs  should  correspond  approximately  to  the  initial 
percentage  of  air  contained  within  the  specimen.  The  initial  volume 
percentage  of  air  contained  in  the  sandy  clay  was  calculated  to  be 
12.7  percent  as  shown  in  Figure  2.3;  the  hydrostatic  pressure  versus 
volumetric  strain  plot  for  Method  V- 5  (Figure  4.1)  shows  a  noticeable 
stiffening  at  a  volumetric  strain  of  approximately  13  percent.  The 
initial  volume  percentage  of  air  for  the  siltstone  specimen  is  2.8 
percent  as  shown  in  Figure  2.4;  although  the  hydrostatic  loading  for 
this  test  was  not  carried  to  a  sufficiently  high  pressure  to  verify 
complete  saturation,  Method  V-5  again  appears  to  indicate  stiffening 
at  the  appropriate  strain,  i.e.  around  2.8  to  2.9  percent.  The  above 
observations,  based  on  air  voids  calculated  for.  two  test  specimens, 
by  no  means  offer  conclusive  evidence  as  to  the  validity  of  Method  V-5; 
but  they  do  illustrate  the  value  of  material-composition  information 
in  the  interpretation  of  mechanical- property  test  data. 

There  is  a  general  trend  of  the  methods  to  combine  into  two  groups : 
one  group,  consisting  of  Methods  V-l,  V-2,  and  V-4,  tends  to  give  low 
strain  values;  and  the  other  group,  consisting  of  Methods  V-3>  V-5, 
and  V-6,  tends  to  give  high  strain  values.  A  possible  bounding  of  the 
volumetric  strain  response  of  soil  under  hydrostatic  loading  may  re¬ 
sult  from  the  use  of  the  standard  Method  V-l  (as  a  lower  bound)  and 
Method  V-3  (as  an  upper  bound) . 

Finally,  without  additional  information  regarding  the  stress  and 
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strain  distribution  within  each  specific  triaxial  specimen,  there 
appears  to  be  no  reason  why  the  standard  Method  V-l,  based  on  actual 
measurements,  should  not  continue  to  be  used  to  develop  plots  for 
data  presentation,  but  the  results  from  this  method,  should  not  be 
used  in  constitutive  property  analyses  without  due  consideration  of 
the  possible  errors  involved. 

4.2  COMPARISON  FOR  SHEAR  TESTS 

Deviator  stress  versus  deviator  strain  plots  for  the  sandy  clay 
specimen  during  the  shear  phase  of  the  test  are  shown  in  Figure  4.3 
for  each  of  the  six  methods  previously  described  for  calculation  of 
deviator  stress  and  strain.  A  table  listing  the  deviator  strains  and 
approximate  shear  moduli  for  two  deviator  stress  ranges  is  included 
in  the  same  figure;  Method  S-l  will  be  considered  the  standard  for 
comparison  purposes.  Note  first  the  similarity  between  all  the  curves 
and  the  relatively  tight  data  band  produced  by  Methods  S-l,  S-2,  S-3, 
and  S-4.  Considering  the  deviator  stress  at  a  deviator  strain  of  4 
percent  as  an  indication  of  ultimate  yield  strength,  the  results  for 
all  six  methods  only  ranged  from  a  high  of  8  percent  greater  (Method 
S-5)  to  a  low  of  7  percent  less  (Method  S- 6)  than  that  given  by  the 
standard  method.  The  shear  moduli  values  calculated  for  the  higher 
deviator  stress  range  (100  to  150  psi)  tend  to  show  more  scatter, 
however.  The  maximum  shear  modulus  in  this  pressure  range  as  calcu¬ 
lated  by  Method  S-5  was  29  percent  greater  them  the  standard  of  2,300 
psi  and  the  lowest  shear  modulus  by  Method  S-6  was  35  percent  less 
than  the  standard. 

The  variation  in  initial  shear  moduli  data  eem  best  be  seen  in 
Figure  4.4,  which  is  an  enlarged  view  of  the  initial  portion  of  the 
deviator  stress-strain  curves  presented  in  Figure  4.3.  Although  there 
is  a  more  noticeable  deviation  from  the  standard  by  Methods  S-5  and 


S-6,  the  tight  banding  of  Methods  S-l,  S-2,  S-3>  and  S-4  is  still  evi¬ 
dent.  The  range  of  shear  moduli  for  the  0-  to  70-psi  pressure  range 
varied  from  a  high  of  115  percent  greater  than  the  standard  for 
Method  S-5  to  a  low  of  48  percent  less  than  the  standard  for  Method 
S-6.  The  deviator  stresses  compared  at  1  percent  deviator  strain 
ranged  from  16  percent  greater  than  the  standard  to  a  low  of  19  per¬ 
cent  less  than  the  standard  value. 

The  deviator  stress-deviator  strain  plot  for  the  shear  phase  of 
the  siltstone  test  is  shown  in  Figure  4.5*  Although  the  curves  tend 
to  produce  the  same  general  trend,  the  values  as  calculated  by  Methods 
S-l,  S-2,  S-3j  and  S-4  are  not  as  close  as  in  the  sandy  clay  example. 
The  deviator  stress  values  taken  at  0.6  percent  deviator  strain  as 
an  indication  of  ultimate  yield  strength  range  from  a  high  of  23  per¬ 
cent  greater  to  a  low  of  21  percent  less  than  the  comparable  deviator 
stress  value  from  the  standard  method.  The  yield  strength  is  not 
shown  in  Figure  4.5,  but  based  on  calculations  for  the  complete  shear 
test  data,  the  maximum  yield  strength  is  approximately  the  same  value 
for  all  methods  although  the  corresponding  deviator  strain  values 
range  from  40  percent  less  than  to  100  percent  greater  than  the  devi¬ 
ator  strain  value  of  1.5  percent  calculated  at  maximum  deviator  stress 
by  the  standard  method.  The  shear  moduli  values  taken  at  high  devia¬ 
tor  stress  levels  (400  to  600  psi)  vary  from  a  high  of  51  percent 
greater  than  the  standard  for  Method  S-5  to  a  low  of  33  percent  less 
than  the  standard  for  Method  S-6.  The  shear  moduli  values  taken  at 
the  medium  deviator  stress  range  (200  to  400  psi)  vary  from  42  percent 
greater  than  to  38  percent  less  than  the  values  calculated  by  the 
standard  method. 

An  expanded  view  cf  the  initial  portion  of  the  deviator  stress- 
deviator  strain  plot  for  the  siltstone  test  is  shown  in  Figure  4.6. 

The  curves  for  this  particular  test,  however,  were  based  primarily 
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on  extrapolation  of  data  measured  at  higher  stress  levels;  and,  al¬ 
though  other  test  data  have  shewn  that  the  various  calculation  methods 
do  affect  the  initial  portion  of  the  shear  results,  use  of  these  par¬ 
ticular  data  will  only  demonstrate  possible  variations.  The  standard 
method  gives  the  minimum  numerical  value  of  deviator  strain  for  any 
given  deviator  stress  level,  while  Method  S-6  still  gives  the  maximum 
numerical  value  of  deviator  strain.  Method  S-5  yields  the  second 
highest  value  of  deviator  strain  up  to  approximately  200-psi  deviator 
stress.  Methods  S-2  and  S-3  give  similar  results.  The  shear  moduli 
values  in  the  0-  to  100-psi  stress  range  decrease  in  value  from  the 
standard  of  555jOOO  psi  to  a  minimum  of  167,000  psi  as  calculated  by 
Method  S-6. 

Several  observations  can  be  made  based  on  the  calculated  stress- 
strain  results  shown  for  the  two  shear  test  examples.  In  general, 
all  the  calculation  methods  will  give  approximately  the  same  values 
for  maximum  yield  strength,  but  significantly  different  values  for 
the  deviator  strain  at  which  maximum  yield  strength  occurs.  Methods 
S-l,  S-2,  S-3,  and  S-4  appear  to  give  approximately  the  same  overall 
stress-strain  results  with  their  curves  bounded  from  above  by  the 
stress-strain  curve  calculated  by  Method  S-5  and  from  below  by  the 
curve  calculated  by  Method  S-6.  Bounding  for  the  initial  shear 
modulus  does  not  appear  feasible  as  yet  since,  at  low  deviator  stress 
values,  any  of  the  methods  could  yield  the  maximum  shear  modulus  ex¬ 
cept  Method  S-6,  which  appears  to  give  the  lowest  values  of  shear 
modulus  regardless  of  material  type.  Method  S-6,  however,  was  based 
on  an  arbitrary  correction  of  axial  deformation  and  the  results  with 
this  method  therefore  can  only  be  considered  as  an  illustration. 

Basically  the  standard  method,  S-l,  is  considered  satisfactory 
at  this  time  for  the  calculation  of  deviator  stress  and  deviator 
strain  from  laboratory  load  and  deformation  measurements  since  it 
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does  tend  to  produce  an  average  of  the  stress-strain  results  as 
calculated  from  the  other  five  methods . 

Although  triaxial  test  results  are  seldom,  if  ever,  presented  in 
terms  of  deviator  stress  versus  volumetric  strain,  comparison  of  such 
plots  based  on  the  various  stress  and  strain  calculation  methods 
previously  described  may  be  beneficial  in  view  of  recent  interest 
by  developers  of  constitutive  models  in  the  volume  change  characteris¬ 
tics  of  soils  and  rocks  during  shear.  Results  for  the  siltstone  ex¬ 
ample  test  are  presented  in  Figure  4.7.  It  is  realized  that  there  is 
some  question  regarding  the  proper  usage  of  Methods  V-5  and  V-6  to 
calculate  volumetric  strain  during  the  shear  phase  of  the  test;  how¬ 
ever,  the  calculated  results  do  illustrate  trends.  The  values  of 
deviator  stress  were  calculated  for  all  six  volumetric  strain  calcu¬ 
lation  cases  by  Method  S-lj  volumetric  strains  were  not  rezeroed  at 
the  start  of  the  shear  phase,  but  were  continued  from  the  last  point 
which  occurred  during  the  hydrostatic  phase. 

Methods  V-l,  V-2,  and  V-4  show  initial  volume  increase  followed 
by  a  volume  decrease  t i  failure.  Method  V-3  shows  an  initial  volume 
increase  followed  by  little  change  in  volumetric  strain  to  peak  devi¬ 
ator  stress.  Method  V-5  shows  little  initial  volume  change  followed 
by  an  increasing  rate  of  decreasing  volume.  Method  V-6  shows  a  rela¬ 
tively  constant  rate  of  increasing  volume  change  followed  by  a  rapid 
increase  near  peak  stress. 

The  principal  observation  that  can  be  made  based  on  the  results 
shown  in  this  example  is  that  any  conclusions  made  from  triaxial  test 
results  as  to  whether  dilation,  compaction,  or  no  volume  change  occurs 
or  does  not  occur  during  shear  of  a  given  material  may  be  completely 
masked  by  the  particular  method  of  data  analysis  chosen  by  the 
experimenter . 
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Figure  4.1  Comparison  of  methods  used  to  calculate 
volumetric  strain  for  the  sandy  clay  test. 
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Figure  k  2  Comparison  of  methods  used  to  calcula 
volumetric  strain  for  the  siltstone  test. 


methods  used  to  calculate  deviator  stress 
le  sandy  clay  test. 
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Figure  4.4  Initial  portion  of  plots  shown  in  Figure  4.3. 
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methods  used  to  calculate  deviator  stress 
le  silt stone  test. 


CHAPTER  5 
CONCLUSION 

The  purpose  of  this  paper  was  to  show  with  specific  examples  the 
possible  influence  exerted  due  to  interpretation  of  raw  load  and  de¬ 
formation  data  on  the  stress-strain  relations  obtained  from  labora¬ 
tory  tests.  The  triaxial  test  was  selected  for  illustration  purposes, 
and  observations  were  made  regarding  physical  response  of  typical 
specimens  subjected  to  both  hydrostatic  loading  and  shear  loading. 

Six  methods  to  calculate  volumetric  strain  and  deviator  stress  and 
strain  were  derived  based  on  a  variety  of  empirical  and  theoretical 
approaches  and  applied  to  calculate  stress- strain  relations  for  two 
widely  different  materials. 

Results  of  the  hydrostatic  phase  calculations  indicated  that  vol¬ 
umetric  strains  were  generally  bound  by  Methods  V-l  (lowest  values) 
and  V-3  (highest  values);  the  actual  hydrostatic  stress-volumetric 
strain  response  of  a  given  soil  or  rock  material  probably  lies  within 
those  bounds.  Results  from  the  shear  phase  calculations  indicated 
that  the  use  of  Method  S-l  to  calculate  the  overall  deviator  stress- 
deviator  strain  curve  is  probably  satisfactory;  however,  selection  of 
an  initial  value  for  shear  modulus  must  depend  to  a  great  degree  on 
the  compatibility  with  moduli  observed  from  tests  conducted  under 
other  states  of  stress. 

Research  is  under  way  to  develop  improved  measurement  systems 
which  attempt  to  circumvent  the  important  influence  of  imposed  bound¬ 
ary  condition  on  the  specimen's  stress-surain  behavior.  Finite  ele- 
ment  analyses  of  laboratory  tests  are  also  under  way  to  provide  some 
insight  into  the  degree  of  influence  of  boundary  conditions  on  the 
distribution  of  stress  and  strain  within  the  specimen  and  possibly  to 
disclose  procedures  for  relating  externally  measured  loads  and 
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deformations  to  the  most  representative  stress-strain  behavior  of  a 
given  test  specimen.  Until  the  results  of  such  research  are  avail¬ 
able,  continued  use  of  the  standard  methods,  i.e.  Method  V-l  for 
volumetric  strain  and  Method  S-l  for  deviator  stress  and  strain,  to 
develop  plates  for  triaxial  test  data  presentation  appears  justified. 

Code  developers  can  provide  very  useful  information  to  experi¬ 
menters  through  their  ability  to  conduct  controlled  parameter  studies 
of  the  influence  of  material  properties  in  ground  shock  calculation; 
results  from  such  studies  can  enable  the  experimenter  to  concentrate 
his  research  efforts  on  those  properties  having  the  greatest  influence 
on  wave  propagation  through  earth  materials.  There  are  certain  favor¬ 
ite  areas  that  experimenters  like  to  investigate  which  may  or  may  not 
have  any  influence  on  the  propagation  of  ground  shock  while  other 
overlooked  areas  might  be  of  prime  importance. 

Hopefully,  through  a  coordinated  research  effort,  such  as  the 
DASA- sponsored  program  at  WES,  significant  results  will  be  attained 
for  use  not  only  in  the  specific  area  of  constitutive  model  formula¬ 
tion  for  ground  shock  codes,  but  also  for  broad  application  in  the 
general  field  of  soil  and  rock  mechanics. 
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